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Abstract 

The  moment  analysis  method  (MA)  has  been  tested  for  the  case  of  2S  — >  2P  ([corejns1  — > 
[corejnp1)  transitions  of  alkali  metal  atoms  (M)  doped  into  cryogenic  rare  gas  (Rg)  matri¬ 
ces  using  theoretically  validated  simulations.  Theoretical/computational  M/Rg  system 
models  are  constructed  with  precisely  defined  parameters  that  closely  mimic  known  M/Rg 
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systems.  Monte  Carlo  (MC)  techniques  are  then  employed  to  generate  simulated  absorp¬ 
tion  and  magnetic  circular  dichroism  (MOD)  spectra  of  the  2S  — >  2P  M/Rg  transition  to 
which  the  MA  method  can  be  applied  with  the  goal  of  seeing  how  effective  the  MA  method 
is  in  re-extracting  the  M/Rg  system  parameters  from  these  known  simulated  systems.  The 
MA  method  is  summarized  in  general,  and  an  assessment  is  made  of  the  use  of  MA  method 
in  the  rigid  shift  approximation  typically  used  to  evaluate  M/Rg  systems.  The  MC-MCD 
simulation  technique  is  summarized  and  validating  evidence  is  presented.  The  simulation 
results  and  the  assumptions  used  in  applying  MA  to  M/Rg  systems  are  evaluated.  The  sim¬ 
ulation  results  on  Na/Ar  demonstrate  that  the  MA  method  does  successfully  re-extract  the 
2P  spin-orbit  coupling  constant  and  Lande  g-factor  values  initially  used  to  build  the  sim¬ 
ulations.  However,  assigning  physical  significance  to  the  cubic  and  non-cubic.  Jahn- Teller 
(JT)  vibrational  mode  parameters  in  cryogenic  M/Rg  systems  is  not  supported. 
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1.  Introduction 

This  study  is  motivated  primarily  by  a  desire  to  understand  and  resolve  the  contradictory  inter¬ 
pretations  of  MCD  spectra  obtained  by  using  the  moment  analysis  method  on  similar  systems  of 
metal  atoms  embedded  in  rare  gas  matrices. 1-3  For  example,  the  work  by  Lund  et  al.1  assigned 
a  spin-orbit  (SO)  splitting  of  -24.5  cm-1  for  the  2P  transition  of  lithium  (Li)  in  solid  argon 
(Li/Ar),  -105  cm-1  for  Li  in  solid  krypton  (Kr),  and  -220  cm-1  for  the  2P  transition  of  Li  in 
solid  xenon  (Xe).  Rose,2,4  working  in  the  same  laboratory,  reported  SO  splittings  of  -327  cm'1 
for  Na/Xe  and  -295  cm”1  for  Li/Xe.  However,  another  paper3  had  contradictory  results,  also 
derived  from  a  moment  analysis,  reporting  an  SO  splitting  of  10  cnr1  for  the  2P  transition  of 
Na/Xe.  Therefore,  a  systematic  investigation  into  the  moment  analysis  method  as  applied  to 
alkali  metal/rare  gas  matrices  (M/Rg)  systems  is  in  order. 

The  aforementioned  papers  use  the  method  of  moments  or  moment  analysis  to  extract  phys¬ 
ical  parameters  from  the  experimental  spectroscopic  data.  These  are  then  used  to  develop 
models  or  theories  of  the  processes  producing  the  values.  If  the  extracted  parameters  are  in¬ 
correct  then  so  too  is  any  interpretation  of  the  data  based  on  those  parameters.  The  moment 
analysis  was  used  in  the  M/Rg  systems  by  analogy  to  f-centers  without  independent  confirma¬ 
tion  of  its  applicability.  This  was  not  necessarily  an  appropriate  analogy.  An  f-center  is  formed 
when  a  negative  ion  is  replaced  by  an  electron  in  a  crystal  lattice;  substantial  interaction  be- 
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tween  the  electron  and  the  surrounding  lattice  is  expected.  However,  the  M/Rg  system  consists 
of  neutral  atoms  in  a  hole  in  a  neutral  lattice;  it  is  not  so  clear  that  a  substantial  interaction, 
as  suggested  by  the  experimental,  moment  analysis  derived,  SO  values,  is  justified.  Therefore, 
the  method  should  be  independently  verified  for  the  system  in  question. 

It  was  found  that  the  theoretically  verified  absorption/MCD  Monte  Carlo  model5  produced 
the  same  basic  absorption  and  MC'D  patterns  observed  experimentally  for  both  spin-orbit  split¬ 
tings  in  the  range  reported  by  the  experimentalists  for  the  Na/Rg  systems  and  using  the  Na  free 
atom  SO  splitting.  Therefore,  it  was  postulated  that  the  moment  analysis  method  as  defined 
for  the  M/Rg  systems  was  biased  to  produce  large  negative  spin-orbit  values.  A  computer  pro¬ 
gram  was  written  to  do  the  full  moment  analysis  evaluation  of  the  simulated  spectra  as  well  as 
tabulated  experimental  spectra.  The  program  was  tested  on  the  published  experimental  spectra 
to  verify  that  it  produced  the  same  values,  within  tabulation  errors,  as  those  published  for  the 
various  parameters.  Once  the  program  was  verified  in  this  way,  it  was  applied  to  the  simu¬ 
lated  spectra.  We  expected  that  both  the  simulated  spectra  produced  using  large  negative  SO 
coupling  and  the  simulated  spectra  using  small  positive  (free-atom  value)  SO  coupling  would 
have  large  negative  moment  analysis  calculated  results  for  the  spin-orbit  coupling,  since  they 
had  the  same  basic  pattern.  However  the  MA  correctly  identified  the  input  SO  coupling  of  each 
simulation. 

2.  Overview  of  the  Method  of  Moments 

The  method  of  moments,  or  moment  analysis  (MA),  provides  a  consistent  way  to  extract  system 
paramters  from  spectra.  Moment  analysis  is  based  on  the  calculation  of  band  moment  integrals 
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over  the  spectral  features,  the  ratios  of  which  are  related  to  physical  quantities.  There  are  two 


types  of  moment  integrals,  initial  and  central.  The  forms  of  the  integrals  are 

Sn  =  J  vnf(v)dv  (2.1) 

for  the  initial  moments  and 

Sn  =  J  (v-  v)n  f  ( v )  dv  (2.2) 

for  the  central  moments  where  v  in  Eqs.  2.1  and  2.2  is  the  photon  frequency,  /  (v)  is  the  spectral 
bandshape  function,  n  is  the  order  of  the  moment,  and  v  is  the  absorption  barycenter.  The 
method  of  moments  can  be  applied  to  any  distribution  function  containing  a  maximum  and 
defined  from  — oo  to  oo.G  As  a  consequence  of  this,  only  discrete  bands  may  be  considered 
using  this  method  since  the  integration  must  take  place  between  points  where  the  band  shape 
function  goes  to  zero.7  In  practice,  this  often  means  that  a  “threshold”  or  background  band 
shape  is  defined  and  subtracted  from  the  spectra  making  it  satisfy  this  requirement.  MA  allows 
different  laboratories  and  research  groups  to  compare  results,  even  if  the  actual  spectra  appear 
quite  different  due  to  equipment  and/or  processing  differences. 

In  the  case  of  IR  (infrared)  spectra,  the  spectrum  is  broken  down  into  sets  of  Gaussian  bands, 
with  integration  over  each  band.6;  8  For  electronic  absorption  and  magnetic  circular  diclrroism 
(MCD)  spectra,  the  integration  is  over  the  entire  band,  with  the  ratios  of  moments  of  various 
orders  from  the  MCD  and  absorption  compared.  An  analysis  similar  to  the  IR  case  may  be 
possible  if  only  the  electronic  absorption  spectrum  is  available. 

Hereafter  all  references  to  moment  analysis  (MA)  refer  to  the  rigid  shift  moment  analysis 
(RSMA)  method.  This  technique  is  widely  used  in  extracting  physical  parameters,  most  notably 
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spin-orbit  coupling,  from  the  MCD  and  absorption  spectra  of  atomic  and  molecular  systems. 
The  method  was  first  applied  to  f-centers,9  but  has  since  been  applied  to  many  diverse  systems, 


including  the  2  S'  — >  2P  absorption  and  MCD  spectra  of  alkali  metals  isolated  in  cryogenic 
rare  gas  solids  (M/Rg).  As  can  be  seen  from  the  moment  equations,  higher  order  moments  are 
increasingly  less  accurate  due  to  the  influence  of  the  wings  of  the  band.'  The  nth  absorption 
and  MCD  moment  integrals  are1 

K  =  j  {^j{S-I)nd£t  (2.3) 

Mn  =  J  (^p)(£-£)nd£,  (2.4) 

where  the  subscript,  n,  indicates  the  order  of  the  moment,  A(£)  is  the  absorption  intensity, 
A  A' (£)  is  the  differential  absorption  of  left  vs  right  circularly  polarized  light,  defined  by 

A A'{£)  =  APCP(£)  -  A'rcp(£),  (2.5) 


where  £  is  the  photon  energy  and  £  is  the  energy  of  the  absorption  barycenter.  The  primed 
terms  indicate  that  the  spectra  are  measured  in  the  presence  of  a  magnetic  field. 


__  fAd£ 
f(f)d£- 


(2.6) 


These  equation  are  of  the  central  moment  type,  with  /  (v)  =  (^)  and  f  (v)  =  f°r 

the  absorption  and  MCD  spectra  respectively.  For  central  moments  such  as  these,  all  odd 
moments  of  completely  symmetric  spectra  are  zero.  Antisymmetric  MCD  spectra  have  zero 
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even  moments. 


The  RSMA  operates  on  the  assumption  that  the  Born-Oppenheimer  (BO)  approximation 
is  obeyed  in  both  the  ground  and  excited  states  and  the  Franck-Condon  (FC)  approximation 
is  sufficient  for  the  electronic  operators  matrix  elements.  For  such  systems,  the  zeroth  and 
first  MCD  moments,  Mo  and  Mi,  contain  all  information  available  from  the  MC'D.  In  systems 
where  the  Jahn- Teller  effect  (JTE)  is  in  operation,  the  BO  approximation  breaks  down  and 
higher  moments  must  be  used.10 

MCD  spectra  are  characterized  by  A,  B,  and  C,  terms.'  In  the  course  of  studying  systems 
in  which  the  MCD  B  term  (See  Table  2.1  for  definition  of  MCD  terms)  is  dominant,  Zgierski 
found  that  the  RSM  equations  generally  used  to  separate  the  A  and  B  term  contributions  of 
the  MCD  spectra  are  inaccurate  when  vibronic  effects  are  present.  The  RSMA  can  be  used 
with  a  vibronically  active  system,'10  but  it  operates  on  the  assumption  that  the  individual 
electronic  transition  moments  between  pairs  of  states  are  vibrational  coordinate  independent 
and  can  be  described  by  a  Franck-Condon  progression.  Zgierski  points  out  that  this  is  not 
a  reasonable  assumption  since  "in  general,  the  adiabatic  potentials  associated  with  different 
states  are  not  parallel  since  electronic  excitation  changes  vibronic  force  fields.  This  implies  that 
the  energy  separation  of  two  coupled  states  is  not  a  constant  but  a  function  of  the  vibrational 
coordinates  (p.  2171)". 11  The  magnetic  mixing  of  vibrations  can  cause  the  MCD  B  term  to 
change  sign  within  the  region  of  a  single,  isolated  allowed  electronic  transition.  Furthermore  the 
first  moment  of  the  MCD  spectrum  is  generally  calculated  with  respect  to  the  center  of  gravity 
of  the  corresponding  absorption  spectrum  and  corresponds  to  the  electronic  A  term.  However, 
when  vibrations  are  present  it  often  is  not  equal  to  the  A  term  "since  the  conventional  definition 
of  the  center  of  gravity  loses  its  applicability  (p.  548)". 12  In  some  cases,  the  magnetic  mixing 
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Table  2.1:  Definitions  of  MCD  Terms 


Term 

Source 

Characteristic  Shape 

A 

Ground  or  exited 
state  degenerate 

+ 

ZA _  1 

A\ 

- 

or  1 

B 

Magnetic  mixing  of 
zero  field  energy  levels 

+ 

+ 

r\ 

or  - 

C 

Ground  state 
degenerate 

+ 

A 

w 

of  vibronic  states  can  simulate  the  presence  of  an  A  term  or  vibronic  coupling  when  only  a  B 
term  is  actually  present.12 

Although  the  2S  — >  2P  transition  of  interest  in  the  case  of  the  M/Rg  systems  contains 
primarily  the  MCD  A  and  C  terms,  consideration  of  the  effects  of  vibrations  and  the  magnetic 
field  on  the  spectra  is  very  important.  It  is  possible  that  matrix  effects  also  contribute  to 
the  mixing  of  the  zero-field  states  making  the  MCD  B  term  and  vibronic  coordinates  more 
important. 

Zgierski  labels  all  deviations  from  the  standard  FC  progression  "non-Condon"  effects.  He 
asserts  that  all  MCD  spectra  are  non-Condon  in  the  sense  that  the  overall  transition  moment 
will  be  coordinate  dependent  if  the  transition  is  induced  by  a  perturbation  that  causes  mixing 
with  a  third  state.  Zgierski  lists  several  references  in  which  he  maintains  it  has  been  shown  that 
non-Condon  effects  are  important  in  vibronically  and  SO  induced  absorption-emission  spectra, 
radiationless  transitions  and  photoelectric  spectra.11 

The  MA  method  is  generally  applied  to  experimental  systems  under  the  implicit  assumption 
that  MA  is  a  valid  way  to  extract  system  properties.  This  study  constitutes  the  first  simulation 
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test / confirmation  of  the  MA  method  for  extracting  the  spin-orbit  (SO)  coupling  of  M /Rg  systems 
[specifically  sodium  (Na)/argon  (Ar)].  The  idea  is  to  build  a  computer  simulated  M/Rg  system 
with  known  SO,  vibronic,  site  symmetry,  and  thermal  parameters  and  known  M-Rg  and  Rg- 
Rg  interaction  potentials  and  then  generate  simulated  2S  — >  2P  electronic  absorption  and 
MOD  spectra  of  this  simulated  M/Rg  system.  These  simulated  spectra  are  then  subjected  to 
a  full  MA  treatment  and  the  SO,  vibronic,  symmetry,  and  thermal  parameters  thus  extracted 
are  compared  to  the  original  parameters  as  an  independent  check  of  the  viability  of  the  MA 
method.  Only  one  other  computer  study  of  a  moment  analysis  method  is  known  to  the  authors. 
R  considers  the  effects  of  spectral  asymmetries  on  the  higher  order  moments.6,8  The  first  of 
these  papers  describes  the  spectral  band  distortions  associated  with  the  detection  equipment  in 
IR  measurements  and  a  technique  for  using  band  moments  to  extract  information.  The  second 
explores  ways  in  which  asymmetries  affect  the  higher  order  moments. 

3.  Specifics  of  MA  for  M/Rg  Systems 

In  the  context  of  this  research  the  MA  method  has  been  applied  to  the  absorption  and  MCD 
spectra  of  metals  trapped  in  cryogenic  rare  gas  solids  [alkali  metal/rare  gas;1-4, 13, 14  alkaline 
earth/rare  gas;15,16  transition  metal/rare  gas17-25].  Those  of  primary  interest  in  high  energy 
density  matter  (HEDM)  research  are  the  alkali  metals,  Li  and  Na.  As  yet  there  has  been 
very  little  experimental  data  on  these  metals  in  the  lighter  rare  gases.  The  Schatz  group  has 
obtained  absorption  and  MCD  spectra  of  Li  in  Ar,  Kr,  and  Xe1  and  Na  in  Xe.2,4  Fajardo  et 
al.  have  shown  that  alkali  metals  can  exist  in  multiple  trapping  sites  in  Rg  matrices  including 
some  novel  tight  trapping  sites  produced  by  the  laser  ablation  method.20  29  Kenney  and  Stowe 
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have  measured  MCD  spectra  of  Li/Ar  tight  trapping  sites  prepared  by  laser  ablation  of  Li.30,31 

After  some  manipulation,  the  results  of  the  different  absorption  and  MCD  moments  can  be 
reduced  to  four  equations  relating  the  zeroth  and  second  absorption  moments,  and  the  first, 
second,  and  third  MCD  moments  to  the  physical  parameters.  They  are1 


M3 

Mi 

A-2 

^0 


-C2  +  3A^  +  -A  2nc 

+  Ac  +  A|rC, 


(3-1) 

(3.2) 

(3.3) 
(3-4) 


where  (  is  the  spin-orbit  parameter  (</  =  (2/3)£,  £  is  the  SO  splitting  of  the  2P  state),  A c  is  the 
cubic  mode  parameter,  A Nq  is  the  non-cubic  mode  parameter,  g  is  the  Lande  (/-factor,  and  k  is 
the  Boltzmann  constant.  Ac  and  A  nc  are  used  to  estimate  the  relative  strengths  of  the  cubic 
and  non-cubic  Jalrn- Teller  (JT)  vibronic  modes  in  the  system.  These  equations  were  originally 
derived  for  f-centers,  therefore  the  inclusion  of  trapping  site  specific/ JT  vibronic  parameters. 

In  practice,  the  spin-orbit  parameter,  (,  is  the  slope  of  the  line  obtained  by  plotting  Mi/Ao 
vs.  — 0.934/2fcT.  The  Lande  (/-factor  is  calculated  from  Eq.  3.1  using  this  (.  The  squared  value 
of  C  obtained  from  the  slope  of  Eq.  3.1  can  then  be  compared  to  that  calculated  from  the  slope 
obtained  by  plotting  M2/Ao  vs.  0.934/4 kT  (Eq.  3.4),  providing  an  independent  check  of  the 
accuracy  of  the  higher  order  MCD  moments. 

For  the  Li/Rg  and  Na/Xe  2S  — >  2P  experimental  spectra1, 2,4  that  have  been  explored 
using  the  MA,  the  values  derived  from  the  MA  equations  for  the  2P  spin-orbit  splitting  were 
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large  and  negative.  The  magnitude  of  the  SO  splittings  obtained  from  MA  for  Li/Ar  and  Li/Kr 
are  substantially  reduced  compared  with  that  of  Li/Xe.1  In  the  Li/Ar  spectra,  the  calculated 
SO  splitting  is  only  a  little  more  than  10%  of  the  Li/Xe  value,  and  the  SO  splitting  in  Li/Kr 
is  approximately  half  the  MA  derived  value  for  Li/Xe.  These  large  negative  values-several 
orders  of  magnitude  larger  and  of  opposite  sign  than  the  free  atom  values-seemed  unphysical. 
The  explanation  of  these  results  was  that  there  was  substantial  interaction  between  the  metal 
atom  valence  electron  and  the  rare  gas.  Although  such  delocalization  of  the  electron  seems 
reasonable  for  the  f-center,  where  an  electron  is  trapped  in  a  negative  ion  vacancy,  it  seems  less 
likely  for  the  case  of  a  metal  atom  in  a  rare  gas.  Furthermore,  simulations  on  the  effects  of 
varying  the  SO  splitting  of  the  impurity  atom  have  shown  that  the  same  gross  spectral  features 
could  be  obtained  using  either  large  negative  spin-orbit  splitting  or  small  positive  spin  orbit 
splitting.  Examination  of  Eq.  3.1  shows  that  for  M-,  /A0  greater  than  zero,  (  is  restricted 
to  either  a  very  small  positive  or  negative  value.  The  experimental  M/Rg  spectra  from  the 
literature  generally  have  roughly  symmetric  absorption,  with  the  exception  of  M/xenon,  and 
MOD  spectra  characterized  by  either  two  down  peaks  (negative)  and  one  up  peak  or  one  down 
peak  and  two  up  peaks  with  increasing  energy.  Consequently,  A0  and  M-\  will  both  have  positive 
signs.  This  tends  to  predispose  the  results  to  require  negative  SO  splitting.  Also,  the  "best 
fit"  to  the  data  described  in  the  literature4  seemed  rather  arbitrary. 

The  observed  spectra  and  MA  results  have  been  attributed1, 2,4  to  the  combination  of  strong 
SO  coupling  and  strong  dynamic  Jahn- Teller  effect  (DJTE)  in  the  2P  manifold.  The  JTE  is  the 
lifting  of  electronic  degeneracies  of  atoms  or  molecules  in  symmetric  trapping  sites  as  the  trapping 
site  is  distorted  by  point-group  specific  vibrations.  The  DJTE  occurs  when  the  trapping 
site  dynamically  accesses  one  or  more  of  the  Jahn- Teller  (JT)  active  vibrations.  Although 
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the  seemingly  too  large  magnitude,  negative  SO  coupling  results  have  been  supported  by  this 
simulation  test  of  the  moment  analysis  method  and  a  simplified  first  principles  calculation,'32 
the  physical  explanations  of  their  source  are  lacking.  The  combination  of  strong  SO  coupling 
and  strong  DJTE  are  contradicted  by  research  by  Ham.33  In  the  absence  of  strong  vibrations, 
SO  coupling  can  be  established  and  a  stable  equilibrium  balancing  the  effects  of  the  neighboring 
Rg  atoms  against  the  orbit  of  the  valence  M  electron  may  be  reached.  When  the  JTE,  or  any 
vibration,  is  present  in  the  solid,  it  becomes  much  more  difficult  to  establish  an  equilibrium 
between  the  two.  Ham  described  the  result  of  this  disruption  of  the  orbit  of  the  electron  as 
"quenching"  the  orbital  angular  momentum.33  The  quenching  of  the  orbital  angular  momentum 
results  in  a  reduction  in  the  magnitude  of  the  SO  coupling.  The  experimental  M/Rg  MA  results 
indicate  an  increase  in  the  magnitude  of  the  SO  coupling,  although  the  sign  has  changed.  If  the 
SO  coupling  is  strong  enough,  it  can  stabilize  the  atom  against  the  JTE,  though  the  magnitude 
of  the  SO  coupling  may  still  be  reduced.33  The  basic  conclusion  that  may  be  drawn  from 
this  is  that  it  is  not  possible  to  have  both  strong  SO  coupling  and  a  strong  JT  interaction. 
Furthermore,  the  JTE  is  defined  as  the  vibronic  lifting  of  electronic  degeneracies  of  atoms  or 
molecules  in  symmetric  sites.  It  is  unlikely  that  all  M/Rg  trapping  sites  can  be  defined  as 
symmetric. 

Another  theory  proposed  to  explain  the  M/Rg  MA  results  invoked  the  external  heavy  atom 
effect  (EHAE)  as  well  as  requiring  symmetric  JT  splitting  of  the  M  J  =  3/2  energy  level  and 
a  12-nearest-neighbor  super  molecule.34  This  theory  has  several  problems.  First  of  all,  the 
EHAE  simply  does  not  apply  to  the  M/Rg  systems.  Although  one  manifestation  of  the  effect 
can  produce  a  "strong  increase  of  the  effective  spin-orbit  coupling,35"  it  tends  to  be  restricted  to 
the  case  of  A  J  =  0,  and  J/0,  a  situation  that  is  not  possible  between  rare  gases  having  integer 
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J  and  alkali  metals  having  half-integer  J.  Furthermore,  the  rare  gases  neon  and  argon  are 
not  heavy  atoms,  yet  the  M/Ar  spectra  are  very  similar  to  those  of  M/Kr  and  M/Xe  systems. 
Na/Ne  absorption  spectra  do  not  exhibit  the  triplet  structure  seen  in  other  M/Rg  combinations, 
however,  it  could  be  that  the  structure  is  simply  not  resolved.27  MCD  spectra  of  Na/Ne  have 
not  been  obtained.  Second,  symmetric  splitting  of  the  M  J  =  3/2  energy  level  is  not  supported 
by  Na/Ar  spectra  and  lastly,  the  12-nearest  neighbor  "super  molecule"  is  very  difficult  to  justify 
in  light  of  the  multiple  2S  — >  2P  absorption  triplets20-29  and  the  various  relative  sizes  of  the 
different  alkali  metal/rare  gas  combinations. 

Most  theoretical  approaches  to  the  MA  question  concerning  the  M/Rg  systems  start  out  from 
a  symmetry  basis:  i.e.  the  site  has  this  symmetry,  so  the  system  will  behave  this  way.  This 
“site”  symmetry  assumption  will  hold  for  F-centers,  or  other  systems  in  which  the  impurity  fits 
comfortably  into  a  single  substitutional  or  interstitial  site  in  a  crystal.  It  may  also  be  useful 
for  studying  atoms  in  molecules,  where  the  surroundings  form  a  specific  and  guaranteed  point 
group.  Site  symmetry  is  not  such  a  good  basis  for  analysis  in  cases  such  as  the  M/Rg  systems 
where  the  M  atom  is  likely  contained  in  a  cavity  larger  than  a  1  Rg-atom  vacancy.  Under  these 
conditions,  identifying  the  site  symmetry,  assuming  that  it  is  consistent  and  symmetric,  is  not 
so  easy.  Furthermore,  it  is  unlikely  that  even  these  assumptions  will  hold  in  all  cases. 

Since  site  symmetry  cannot  be  guaranteed  in  many  systems,  including  most  M/Rg  systems, 
the  question  “must  the  site  be  symmetric  for  MA  to  be  valid”  becomes  significant.  It  seems, 
based  on  careful  reading  of  the  original  moment  analysis  papers,9,30  the  papers  on  the  physical 
basis  for  the  model,3 1 ' 38  and  our  experience  with  the  simulations,  that  the  answer  to  this  question 
is  “NO.” 

In  fact  our  simulations,  described  below,  have  consistently  and  conclusively  demonstrated 
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that  site  symmetry  does  not  play  a  significant  role  in  the  MA  derived  spin-orbit  results.  The 


importance  of  this  lack  of  variation  with  trapping  site  symmetry  in  the  MA  spin-orbit  and  g- 
factor  results  is  that  it  shows  there  is  no  reason  for  any  assumption  of  site  symmetry  in  evaluating 
a  system  using  the  MA. 

4.  Monte  Carlo  simulation  of  the  absorption  and  MCD  spectra  of  M/Rg 

The  simulation  process  utilized  in  modeling  the  2S  — >  2P  transition  in  a  M/Rg  system  uses 
a  classical  Monte-Carlo  (MC)  scheme  based  upon  the  original  Metropolis  et  al.  algorithm39 
exactly  as  implemented  by  Boatz  and  Fajardo  for  the  case  of  electronic  absorption  spectroscopy40 
with  extension  to  include  MCD  spectroscopy  accomplished  by  Kenney  et  al.5  A  MC-MCD 
simulation  begins  by  choosing  an  initial  M-Rg  configuration  {M,  R1;  R2, ...,  Rjv /initial,  typically, 
representing  an  idealized  trapping  site  [e.g.,  M  in  a  single  substitutional  site  (one-atom  vacancy) 
with  Oh  symmetry].  The  full  MCD  spectrum  for  a  2S  — >  2P  transition  in  a  M/Rg  system  in 
a  specific  M-Rg  configuration  {M,  Rj,  R2, ...,  Rn}  is  a  "stick  MCD  spectrum"  consisting  of  12 
lines,  appropriately  placed  on  the  x  (energy)  axis  at  transition  energies  hvfi,  whose  magnitudes 
and  directions  (+  or  -)  on  the  y  (AA7)  axis  are  determined  by  computing  Eq.  2.5  as  restated 
here, 

AA7  ( Ufi )  =  A'lcp  (vfi)  —  A'rcp  (vfi) ,  (57) 

for  all  possible  choices  of  i  and  /  and  the  ( 7  )  indicates  the  measurement  is  taken  in  a  magnetic 
field.  To  simulate  a  real  MCD  spectrum,  many  stick  MCD  spectra,  each  arising  from  a  different 
Rg  configuration,  are  appropriately  averaged  together.  The  MC  energy  optimization  scheme 
operates  in  the  2S  ground  state  manifold.  At  this  stage,  the  very  small  Zeeman  perturbation 
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of  the  2S  manifold  is  neglected.  MCD  stick  spectra  arising  from  more  favorable  configurations 
are  weighted  more  heavily  in  the  averaging  process  than  those  arising  from  the  less  favorable 
configurations.  In  actual  practice,  the  LCP  and  RCP  contributions  to  the  MCD  spectrum  ap¬ 


pearing  in  Eq.  2.5  are  accumulated  and  stored  separately.  This  allows  the  electronic  absorption 
spectrum  to  be  recovered  from  the  MC-MC'D  simulation  as' 


A'  —  A  —  ~[A'lcp  +  A'rcp\. 


(4.1) 


Eq.  4.1  expresses  the  fact,  well  established  experimentally,  that  the  electronic  absorption 
spectra  of  M/Rg  systems  are  essentially  unaffected  by  the  application  of  an  external  magnetic 
field. 


4.1.  Orientational  averaging  in  MC-MCD  simulations: 

In  an  MCD  experiment,  the  propagation  direction  of  the  light,  which  is  always  parallel  to  the 
B  vector,  rigorously  defines  the  laboratory  frame  z  axis.  The  laboratory  frame  x  and  y  axes 
are  set  parallel  to  the  sapphire  deposition  window  and  perpendicular  to  B  with  the  origin  at 
the  M  nucleus.  All  Rg  positions,  eigenvectors,  and  selection  rules  are  described  with  respect  to 
this  laboratory  frame.  To  model  the  effects  of  arbitrary  Rg  lattice  orientation  relative  to  the 
laboratory  frame,  the  2P  interaction  matrix5  is  separated  into  M/Rg  and  spin-orbit  +  Zeeman 
matrices 

[M/Rg  +  SO  +  Zeeman]  =  [M/Rg]  +  [SO  +  Zeeman].  (4.2) 

Arbitrary  orientations  of  the  Rg  matrix  relative  to  the  laboratory  frame  can  be  modeled  by 
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re-writing  Eq.  4.2  as 


[M/Rg  +  SO  +  Zeeman]  =  [J7][M/Rg][£/]  1  +  [SO  +  Zeeman],  (4.3) 

where  U  is  an  arbitrary  randomly  generated  unitary  rotation  matrix  (e.g.,  an  Eulerian  rotation 
matrix).  The  MC-MCD  simulation  process  is  carried  out  as  before  using  Eq.  4.3  rather  than 
Eq.  4.2. 

4.2.  Lattice  temperature  vs.  magnetic  temperature  in  MC  MCD  simulations: 

Temperature  effects  appear  in  two  distinct  places  in  MC-MCD  simulations.  A  lattice  temper¬ 
ature,  T,  must  be  specified  for  the  MC  simulation  process.  It  may  be  argued  that  a  proper 
representation  of  quantum  mechanical  zero  point  vibrational  motions  of  the  M/Rg  lattice  at 
temperature  T  requires  the  classical  MC  simulation  to  be  carried  out  at  a  classical  lattice  tem¬ 
perature  T'  >  T.41  However,  a  temperature  for  the  Boltzmann  factor  of  the  Zeeman-split  2S 
term  also  must  be  included.  This  "magnetic  temperature"  always  should  be  the  actual  system 
temperature,  T.  The  MC-MCD  simulation  code  has  the  option  of  specifying  separate  values  for 
the  lattice  temperature  T'  and  the  magnetic  temperature  Tmag .  Lawrence  and  Apkarian  discuss 
the  relationship  between  T'  and  Tmag  in  classical  MC  simulations  of  doped  Rg  systems.41 

4.3.  Truth-testing  of  MC-MCD  Simulations: 

At  this  point  it  should  be  noted  that  the  MC-MCD  simulation  code  based  upon  the  theory 
outlined  in5  and  used  with  this  work  has  passed  the  following  tests  in  Na/Ar.  The  MC-MCD 
spectra  (and  associated  absorption  spectra)  (1)  exhibit  MCD  and  absorption  band  profiles  that 
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look  very  similar  to  the  experimental  spectra  of  M/Rg  systems,  (2)  the  MCD  and  absorption 
band  profiles  properly  overlay  each  other  in  energy,  (3)  properly  invert  when  the  direction  of 
the  magnetic  field  is  reversed,  (4)  increase  or  decrease  in  amplitude  with  increases  or  decreases 
in  B,  (5)  yield  proper  spin-orbit  coupling  constant  magnitudes  and  signs  from  band  moment7 
analyses,  (6)  are  suitably  responsive  to  temperature  effects  and  orientational  averaging. 

5.  Evaluation  of  the  Simulations  Using  the  MA  Method 

A  computer  program  was  written  to  calculate  the  moment  integrals  on  the  simulated  spectra 
produced  in  the  MC-MCD  program.  Since  the  2P  SO  coupling  constant  sign  and  magnitude, 
trapping  site  symmetry,  magnetic  field,  and  temperature(s)can  be  chosen  explicitly  in  the  sim¬ 
ulation,  a  successful  MA  method  must  recover  these  parameters  from  the  simulated  2 S  — >  2P 
absorption  and  MCD  spectra  that  are  calculated  to  arise  from  this  simulated  M/Rg  system. 
Therefore,  if  the  MA  method  did  not  produce  the  expected  spin-orbit  splitting,  it  would  be 
precisely  falsified  for  the  M/Rg  systems.  The  MA  consistently  reproduced  the  input  spin-orbit 
splitting  to  within  reasonable  variations.  As  yet,  the  only  times  in  which  the  MA  has  not 
reproduced  the  input  SO  splitting  within  the  error  bounds  is  when  the  input  magnetic  held  is 
large.  In  the  cases  examined,  as  the  magnetic  held  is  increased,  the  MA  increasingly  underesti¬ 
mates  the  SO  splitting.  This  is  not,  however,  a  major  concern,  since  most  experimental  MCD 
spectra  are  acquired  with  magnetic  helds  much  smaller  than  those  produced  these  deviations  in 
the  simulations. 

In  this  study,  a  MA  data  set  is  constructed  by  generating  a  set  of  simulated  absorption  and 
MCD  spectra  using  multiple  magnetic  temperatures;  typically,  Tmag  =  4,  6,  8,  10,  20,  and  30 
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K .  are  used.  A  single  substitutional  site  was  used  in  generating  all  of  the  simulated  absorption 
and  MCD  spectra.  The  moment  integrals,  Eqs.  2.3  and  2.4,  are  evaluated  numerically  using 
Simpson’s  rule.  For  each  temperature,  the  first,  second,  and  third  MCD  moments  and  the 
zeroth  and  second  absorption  moments  are  calculated.  The  £2  term  used  in  calculating  the 
cubic  and  non-cubic  vibronic  mode  parameters,  and  A2NC,  is  derived  from  the  slope  of  the 
spectral  data,  though  it  can  also  be  calculated  directly  from  (.  Then  g,  (.  A 2C  and  A2NC  are 
calculated  from  Eqs.  3.1,  3.2,  3.3,  and  3.4.  The  technique  was  tested  by  tabulating  spectra 
from  the  literature  and  comparing  the  results  with  those  published.  The  SO  splittings  obtained 
from  the  MA  generally  were  in  good  agreement  with  the  values  input  into  the  simulations,  with 
deviations  accounted  for  by  the  numerical  errors  associated  with  tabulating  the  spectra  and 
calculating  the  moment  integrals. 

Single  spectrum  MA  of  simulation  or  experimental  data  is  also  possible.  This  requires  that 
an  assumption  for  the  g-factor  be  made  for  the  initial  calculation.  For  simulation  data  a  g- 
factor  of  1,  the  average  of  the  atomic  P-state  g- factors,  is  a  good  approximation.  This  is  because 
unaltered  simulation  spectra  consistently  produce  g  «  1.00  in  the  MA  evaluation. 

For  experimental  data,  the  choice  of  g  is  not  so  clear,  since  it  is  conceivable  that  the  true 
value  may  differ  from  the  atomic  average.  In  fact,  Adrian4"  claims  that  the  matrix  g-value 
is  necessarily  smaller  than  the  free  atom  g-value  due  to  mixing  of  the  matrix  atom  properties 
into  the  isolated  atom  wave  function.  This  overlap  causes  the  matrix  SO  interaction  to  have 
the  effect  of  reducing  the  g-value,  regardless  of  the  nature  of  the  trapping  site,  although  the 
difference  is  ~10~4  or  less  for  a  trapped  hydrogen  (H)  atom.42  Reductions  in  the  magnitude 
of  the  g-value  were  also  calculated  by  Smith  [in,43  p.  25]  for  the  alkali  atoms  trapped  in  solid 
rare  gases. 
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Therefore,  when  the  experimental  spectrum  tested  is  from  published  sources,  the  published 
5-value  may  be  used.  If  no  5-factor  is  available,  or  its  accuracy  may  be  in  doubt,  an  estimate 
of  the  5-factor  can  usually  be  obtained  from  the  ratios  of  Mi/Aq  and  M^/Aq.  M2/A0  gives 
an  estimate  of  £2,  from  which  the  positive  and  negative  roots  may  be  substituted  into  Eq.  3.1, 
producing  two  possible  5-values.  Generally  one  of  these  can  be  rejected,  identifying  the  correct 
root  of  <j2  and  therefore  the  best  estimate  of  the  SO  splitting,  £.  This  technique  for  estimating 
5  fails  when  ( M2/A0 )  <  0.  This  unphysical  situation  has  often  occurred  in  simulation  results. 
It  is  not  known  if  it  may  also  happen  with  experimental  results,  though  some  researchers  have 
not  used  this  value  in  analyzing  their  data  since  it  is  unreliable.1 

Since  the  single  spectrum  MA  technique  described  here  rests,  to  a  certain  extent,  on  the 
choice  of  5- factor,  it  seemed  prudent  to  test  the  sensitivity  of  the  calculation  to  it.  It  was 
found  that  the  5-factor  has  no  effect  on  the  cubic  mode  parameter,  A^,  and  only  limited  effect, 
through  the  value  of  <j2  used,  on  the  non-cubic  mode  parameter,  A?NC.  The  same  cannot  be 
said  of  the  SO  splitting.  Our  simulation  tests  indicate  that  the  choice  of  5-value  affects  the  MA 
calculated  SO  splitting  in  a  consistent  pattern,  with  the  calculated  SO  splitting  increasing  with 
increasing  5-value.  For  systems  having  large  SO  splittings,  the  variation  in  the  SO  splitting 
with  5-value  is  significant,  but  not  extreme.  When  the  SO  splitting  is  small,  however,  these 
variations  can  be  greater  than  the  total  SO  splitting. 

In  this  test,  the  5- value  used  to  calculate  the  SO  splitting  was  varied  from  0.0  to  2.0  in  steps 
of  0.5.  In  each  case,  the  calculated  SO  splitting  was  reduced  by  ~  20  cm-1  for  5  =  0.0,  and 
increased  by  ~  20  cm-1  for  g  =  2.0  compared  to  the  value  calculated  for  5  =  1.0,  the  free  atom 
value.  The  trapping  site  was  a  one-atom  vacancy,  the  magnetic  field,  B  =  IT,  and  the  selected 
magnetic  temperature,  Tmag  =  10  K.  Table  5.1  displays  the  results  of  this  test.  The  spin  orbit 

19 


This  page  is  Distribution  A:  Approved  for  public  release,  distribution  unlimited 


Table  5.1:  Effect  of  the  g- value  on  the  single  spectrum  MA  calculated  SO  coupling. 


g 

Calculated  SO  splitting 

-170.0  cm  1 

170.0  cm  1 

17.0  cm  1 

0.0 

149.1 

-3.9 

im 

159.6 

6.5 

1.0 

-169.9 

170.0 

16.9 

1.5 

-159.4 

180.4 

27.4 

2.0 

-149.0 

190.8 

37.8 

splittings  above  the  double  line  are  the  input  values  used  in  generating  the  simulated  spectra. 
The  values  listed  across  from  the  g-values  are  those  calculated  using  that  g-value  and  the  ratio 


Mi/A0. 


6.  Exploration  of  Governing  Factors  in  MA  Results 

6.1.  Spin-orbit  splitting  and  the  MA  method 

In  general,  the  absorption  spectra  for  all  simulation  tested  spin-orbit  coupling  values  are  very 
similar.  In  fact,  when  the  SO  splitting  is  in  the  range  ±  50  cm-1,  the  same  trapping  site 
at  the  same  temperature  produces  identical  absorption  spectra  down  to  expected  statistical 
errors  inherent  in  Monte  Carlo  (MC)  sampling.  The  exceptions  tested,  ±170  cm”1,  have  SO 
splittings  on  the  same  order  of  magnitude  as  the  M-Rg  and  Rg-Rg  matrix  interaction  energies. 
The  diatomic  interaction  potentials  show  that  the  argon-argon  binding  energy  is  around  100 
cm-1,  and  the  sodium-argon  (Na-Ar)  binding  energy  is  around  58  cm-1. 

This  suggests  that  the  MA  results  depend  primarily  on  the  MCD  spectral  features.  To 
explore  this  theory,  several  tests  were  instituted.  In  the  first  test,  two  sample  data  sets  were 
constructed.  One  of  these  had  the  absorption  spectrum  generated  with  input  SO  splitting  at 
-170  cm-1  and  the  MCD  spectrum  generated  with  SO  splitting  of  17  cm”1;  the  second  had  the 
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reverse  with  the  absorption  SO  splitting  =  17  cm-1  and  MC'D  SO  splitting  =  -170  cm”1.  These 
input  SO  splitting  values  were  chosen  since  both  produce  an  MOD  spectrum  characterized  by 
two  down  peaks  and  one  up  peak  with  increasing  energy.  This  test  strongly  supported  the 
preeminence  of  the  MOD  spectrum  in  determining  the  MA  SO  splitting.  For  the  set  with  the 
MOD  generated  using  an  input  SO  splitting  of  17  cm”1  the  MA  calculated  SO  splitting  was 
16.87  cm”1,  easily  in  the  range  of  expected  MA  SO  splitting  values  for  a  simulation  set  with 
input  SO  splitting  of  17  cm”1.  In  the  reverse  test  where  the  MOD  SO  splitting  was  -170  cm”1, 
the  MA  calculated  SO  splitting  was  -263.43  cm”1.  The  reason  for  the  substantial  increase 
in  the  magnitude  of  the  SO  splitting  in  this  case  is  not  entirely  understood,  but  is  probably  a 
function  of  the  relatively  large  contribution  of  the  outer  edges  of  the  MOD  spectra  in  comparison 
to  the  absorption,  as  supported  by  the  tests  of  the  spectral  bandwidth  reported  later  in  this 
paper.  Typically,  the  relative  intensity  of  the  MOD  spectrum  is  several  orders  of  magnitude  less 
than  the  absorption  spectrum  over  the  entire  spectral  range.  In  this  case  however,  the  MCD 
band  was  broader  than  the  absorption  band,  consequently  it  has  a  significant  intensity  where 
the  absorption  is  small.  This  could  occur  experimentally  if  the  zero  point  of  the  absorption  is 
incorrectly  defined  or  absorption  is  present,  but  not  significant  compared  to  the  maximum  of  the 
spectra.  The  simulations  have  shown  that  A  A!  I A  is  larger  for  larger  magnitude  SO  splittings 
than  for  smaller  magnitude  SO  splittings.  Therefore  if  A  A I A  is  artificially  inflated,  the  result 
is  a  larger  magnitude  calculated  SO  splitting,  with  the  peak  pattern  determining  a  negative 
sign  for  the  SO  splitting.  The  magnitude  of  A  A I A  also  plays  a  role  in  whether  the  MCD 
peak  pattern  of  down,  down,  up  is  calculated  as  positive  SO  coupling  or  negative  SO  coupling. 
Smaller  A  A I A  tend  to  produce  positive  SO  splitting  values  for  this  pattern.32 

To  further  explore  the  effect  of  the  MCD  on  the  MA  calculated  SO  results,  a  second  test 
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in  which  the  MCD  of  a  specific  set  is  multiplied  by  a  scale  factor  was  conducted.  In  this  test 
the  resultant  MA  calculated  SO  splitting  was  consistently  equal  to  the  value  of  the  input  SO 
splitting  multiplied  by  the  scale  factor.  Similarly  the  calculated  g-factor  varied  directly  with 
the  scale  factor.  When  no  scaling  occurred,  the  MA  calculated  SO  splitting  is  approximately 
equal  to  the  input  SO  splitting,  and  the  g-factor  is  approximately  equal  to  the  averaged  atomic 
M  P-state  value  of  1.  Although  it  may  be  very  difficult  to  experimentally  verify,  this  suggests 
that  faulty  MA  results  can  occur  if  something  other  than  the  SO  splitting  is  a  major  factor 
in  producing  the  relative  heights  of  the  MCD  and  absorption  bands.  It  further  suggests  that 
the  g-factor  might  provide  an  estimate  of  the  degree  of  error  in  the  SO  splitting  measurement. 
Table  6.1  demonstrates  the  variation  in  the  MA  calculated  SO  splitting  and  calculated  g-factor 
when  the  MCD  spectra  for  a  1-atom  vacancy  with  input  SO  splitting  of  17  cnr1  is  scaled  by 
multiplying  the  MCD  intensity  by  a  selected  scale  factor.  Included  in  the  table  for  comparison 
are  the  input  and  MA  calculated  SO  splittings,  and  the  MA  calculated  g-factors.  Notice  that 
the  g-factors  for  the  unsealed  sets  (column  6,  Table  6.1)  are  all  approximately  1.  Table  6.2  is 
the  corresponding  table  for  input  SO  coupling  of  -170  cm-1.  The  values  here  can  be  compared 
with  those  in  columns  3-6  in  Table  6.1. 

The  possibility  that  some  factor  might  influence  the  width  of  the  MCD  bands  with  respect  to 
the  absorption  bands  was  also  considered.  The  spectral  bandwidths  of  the  simulated  absorption 
and  MCD  spectra  were  separately  broadened  or  narrowed  by  adjusting  energy  steps  (a.k.a. 
nominal  simulation  bin  size)  in  the  calculation  of  the  moment  integrals.  In  some  cases,  the 
energy  barycenter  of  the  absorption  or  MCD  was  shifted  with  respect  to  the  original  placement. 
It  was  found  that  broadening  or  narrowing  the  absorption  spectra  has  very  little  effect  on  the 
SO  splitting  results  (see  Table  6.3.)  However,  the  same  is  not  true  of  changes  in  the  MCD  band 
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Table  6.1:  Effect  of  scale  factor  on  spectra,  initial  SO  splitting  =  17  cm 


17  cm  1 

from  scaled 

from  original 

Scale  factor 

Calc.  SO 

g-factor 

Input  SO 

Calc.  SO 

g-factor 

1.0 

16.87 

1.00 

17  cm-1 

-1.0 

-16.87 

-1.00 

-17  cm-1 

-16.89 

1.00 

10.0 

168.72 

9.97 

170  cnr1 

167.74 

0.96 

-10.0 

-168.72 

-9.97 

-170  cm-1 

-169.02 

1.07 

5.0 

84.36 

4.99 

N/A 

N/A 

-5.0 

-84.36 

-4.99 

N/A 

N/A 

3.0 

50.61 

2.99 

50  cnr1 

49.66 

0.99 

-3.0 

-50.61 

-2.99 

-50  cm-1 

-49.68 

1.00 

0.6 

10.12 

0.60 

10  cm-1 

9.94 

1.00 

-0.6 

-10.12 

-0.60 

N/A 

N/A 

Table  6.2:  Effect  of  scale  factor  on  spectra,  initial  SO  splitting  =  -170  cm  1. 


170  cm  1 

from  scaled 

Scale  factor 

Calc.  SO 

g-factor 

0.1 

-16.80 

0.10 

-0.1 

16.80 

-0.10 

0.3 

-50.40 

0.31 

-0.3 

50.40 

-0.31 

1.0 

-168.00 

1.03 

-1.0 

168.00 

-1.03 

width.  These  had  a  proportionate  effect  on  the  MA  calculated  SO  splitting.  Broadening  or 
narrowing  either  the  absorption  or  MOD  spectrum,  but  not  the  other,  in  all  cases  significantly 
affected  the  cubic  and  non-cubic  vibronic  mode  parameters,  sometimes  making  one  or  more 
of  the  M2/A0  ratios  less  than  zero,  an  unphysical  result.  Broadening  or  narrowing  both  the 
absorption  and  MOD  spectra  tended  to  keep  these  values  more  stable.  Table  6.4  gives  the 
energy  bin  sizes  and  starting  energies  for  the  MOD  and  absorption  spectra.  The  unaltered 
(default)  bin  size  is  10  cm-1. 

Figure  6.1  shows  the  effect  of  varying  the  MC'D  bin  size,  plotting  the  SO  splitting  from  Table 
6.3,  against  the  variation  of  the  MCD  bin  size  with  respect  to  the  simulation  bin  size.  Although 
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Table  6.3:  The  effects  of  spectral  bandwidth  on  MA  results.  Tests  marked  with  an  ‘*’  have  the 
MC'D  or  absoprtion  energy  center  shifted  with  respect  to  the  original  placement.  The  Lande 
g- factor  calculated  from  Aa.  is  “NAN"  (not  a  number)  when  ^  is  less  than  zero  for  at  least 
one  simulation  in  the  data  set. _ 


SO  coupling 

Lande  g-factor 

Squares  of  (Vibronic) 

(from  slope) 

Calculated  from: 

Cubic  Mode 

Non-Cubic  Mode 

Test 

(cm”1) 

slope 

M2 

aq 

input  SO 

Param.  (A;?) 

Param.  (A^c) 

1 

16.87 

0.999 

395.525 

71991.039 

2 

18.56 

1.097 

NAN 

22297.029 

43489.871 

3 

15.18 

— 

0.965 

-21336.740 

101785.312 

4 

16.87 

EMM 

1.038 

-23627.068 

120000.828 

5 

16.87 

0.913 

19888.279 

32861.023 

6 

20.26 

1.197 

NAN 

44570.574 

15713.099 

7 

13.49 

na 

0.863 

-43943.586 

134465.422 

8* 

15.08 

EH 

NAN 

-20047.752 

100496.984 

9* 

15.18 

0.899 

345.754 

58283.531 

10* 

16.98 

1.003 

NAN 

30884.021 

21864.361 

11 

16.86 

1.059 

-52429.555 

177575.953 

12 

25.31 

1.496 

1.495 

mmm 

561.952 

162361.500 

13 

25.25 

1.493 

NAN 

mstwm 

112872.141 

-64713.969 

14* 

17.73 

1.052 

NAN 

mmm 

295392.812 

-229599.203 

15* 

19.96 

1.174 

NAN 

1.002 

793403.500 

-727628.562 

there  are  minor  variations  in  the  calculated  SO  splitting  depending  on  the  difference  between 
the  entered  values  for  the  MC'D  and  absorption  bin  sizes  or  whether  or  not  the  starting  value 
of  the  MC'D  and  absorption  energies  have  been  adjusted  to  keep  the  energy  barycenter  in  its 
original  position,  it  is  safe  to  say  that  the  MOD  bin  size  is  the  dominant  factor  in  the  calculated 
SO  splitting. 

Interestingly,  small  displacements  of  the  MOD  with  respect  to  the  absorption  do  not  have  a 
significant  effect  on  the  MA  calculated  SO  splitting.  Sample  data  files  were  created  by  shifting 
the  MC'D  and  absorption  spectrum  energy  axes  with  respect  to  one  another  by  a  small  interval 
5e,  then  MA  was  done  on  a  “data”  set.  The  unshifted  data  set  had  a  MA  calculated  SO  coupling 
of  16.87  cm”1.  Shifting  the  MC'D  relative  to  the  absorption  100  cm”1  made  the  MA  calculated 
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Table  6.4:  Spectra  broadening  test  parameters.  Tests  marked  with  used  the  same  bin  size 
for  the  absorption  and  MCD. 


TEST 

MCD 

Absorption 

Energy  bin 
size  (cm-1) 

Starting  Energy 
(cm-1) 

Energy  bin 
size  (cm-1) 

Starting  Energy 
(cm-1) 

1 

10 

13990 

10 

13990 

2 

11 

13594 

10 

13990 

3 

9 

14386 

10 

13990 

4 

10 

13990 

11 

13594 

5 

10 

13990 

9 

14386 

6 

12 

13198 

10 

13990 

7 

8 

14782 

10 

13990 

8 

9 

14440 

10 

13990 

9* 

9 

14440 

9 

14440 

10 

10 

13990 

9 

14440 

11 

10 

13990 

12 

12* 

15 

12040 

15 

13 

15 

12040 

10 

13990 

14 

11 

10 

13990 

15 

11 

10 

13990 

SO  coupling  17.05  cnr1  and  16.68  cm-1,  when  the  MCD  was  shifted  to  a  lower  energy  or  a 
higher  energy,  respectively,  compared  to  the  absorption.  Considering  the  magnitude  of  the  total 
relative  shift  between  the  MCD  and  absorption,  these  changes  in  the  calculated  SO  splitting  are 
negligible. 

Non-spin-orbit  sources  that  may  contribute  to  the  MCD  spectra  include  the  pseudo-Jahn- 
Teller-effect  (PJTE)  and  orbit  lattice  coupling.  The  PJTE  is  the  vibrational  mixing  of  "the 
finite  electronic  state  with  another  closely  located  electronic  level",44  which  can  influence  the 
energies  and  spacing  of  the  levels.  Orbit-lattice  couplings,  coupling  of  the  electron  orbit  to 
lattice  vibrations,  may  sizably  change  the  energy  spacings  of  the  levels.45  It  has  also  been 
suggested  that  external  axial  crystal  fields  can  mimic  SO  coupling.46 

The  basic  conclusion  that  can  be  drawn  from  the  tests  of  the  MA  governing  factors  is  that 


25 


This  page  is  Distribution  A:  Approved  for  public  release,  distribution  unlimited 


Calculated  SO  splitting  (cm1) 


Figure  6.1:  The  effect  of  adjusting  the  MCD  bin  size  on  the  calculated  spin  orbit  splitting 
parameter.  <,  >,  and  =  in  the  legend  means  the  MCD  bin  size  is  set  smaller  than,  larger  than, 
equal  to,  respectively,  the  corresponding  absorption  bin  size;  c  means  the  MCD  bins  are  centered 
on  the  true  simulation  absorption  barycenter,  n  means  they  are  not. 

while  the  MA  method  is  relatively  impervious  to  any  changes  to  the  absorption  spectra,  the 
same  cannot  be  said  of  variations  in  the  MCD  spectra.  Depending  on  the  source  of  variation 
the  effects  of  modifications  to  the  MCD  spectra  can  range  from  negligible  to  substantial. 

6.2.  The  Jahn- Teller  vibronic  mode  parameters,  and  A^c 

While  the  MA  has  been  extremely  successful  in  extracting  the  spin-orbit  splitting  of  the  simu¬ 
lations,  and  hence  the  g-factor,  the  results  for  the  square  of  the  cubic  vibronic  mode  parameter 
and,  to  a  lesser  extent,  the  square  of  the  non-cubic  vibronic  mode  parameter,  have  been  highly 
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erratic.  Dramatic  changes  in  the  value  of  the  square  of  the  cubic  mode  parameter  occur  with 
simple  changes  in  the  choices  of  magnetic  temperatures  evaluated  in  a  single  set  of  simulations. 
In  some  cases,  this  value  has  been  negative,  which  is  not  physically  possible,  since  it  is  a  squared 
parameter.  Since  these  parameters  are  vibronic  in  nature,  and  specifically  apply  to  symmet¬ 
ric  sites,  their  inclusion  in  MA  of  M/Rg  systems,  with  the  possible  exception  of  Li/Xe,  are 
suspect.  They  are  a  holdover  from  the  f-center  origins  of  the  MA  equations.  Consequently 
their  applicability  to  the  M/Rg  systems  was  not  questioned  or  tested.  Although  the  failure  of 
vibronic  modes,  Aq  and  A2NC,  in  simulation  results  is  not  proof  that  the  values  are  not  valid 
in  experimental  results,  it  certainly  casts  doubt  on  their  applicability.  Furthermore,  dynamic 
Jahn- Teller  (DJT)  analysis  of  Li/Xe  results4'  also  cast  doubt  on  the  validity  of  the  MA  derived 
values  of  the  vibronic  modes,  A/,  and  A2NC.  In  this  study,  a  Jahn- Teller  analysis,  assuming 
an  octahedral,  Oh,  trapping  site  and  equal  coupling  of  the  t2g  and  eg  vibrations,  was  able  to 
replicate  the  observed  absorption  and  MC-D  spectra  with  reasonable  success.  However,  the 
values  of  and  A2NC  entered  into  the  JT  analysis  producing  these  results  were  NOT  the 
values  obtained  through  the  MA.  Rather  the  parameters  were  adjusted  until  the  best  fit  to 
the  data  was  obtained.  Furthermore,  the  theory  is  insufficient  to  account  for  other  alkali  atom 
spectra.47 

Figure  6.2  is  a  graphical  representation  of  the  scatter  in  the  cubic  mode  parameter  as  it  varied 
with  the  method  test  shown  in  Table  6.3.  The  values  from  tests  14  and  15  were  excluded  from 
this  table  as  they  were  so  large  that  they  made  the  range  scale  unreasonable.  Although  there 
is  an  apparent  trend  to  higher  values  of  the  cubic  mode  parameter  with  increasing  apparent  SO 
splitting,  this  cannot  be  assumed  to  be  the  case;  by  focusing  only  on  the  values  right  around  17 
cm-1,  it  can  be  seen  that  the  distribution  is  very  random.  Negative  values  of  the  cubic  mode 
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parameter  are  not  possible  since  it  is  a  squared  value,  therefore  all  values  below  the  blue  line 
represent  unphysical  results.  The  non-cubic  mode  parameter  is  not  plotted,  but  typically  scales 
inversely  with  the  cubic  mode  parameter,  going  negative  (again  an  unphysical  result)  when  the 
cubic  mode  parameter  is  very  large.  In  all  of  these  cases,  the  same  trapping  site,  a  single 
substitutional  site  (a.k.a.  one-atonr  vacancy),  was  used. 

The  Jahn- Teller  effect  is  defined  to  be  the  vibronic  lifting  of  electronic  degeneracy  in  sym¬ 
metric  trapping  sites.  If  the  trapping  site  is  inherently  asymmetric,  then  this  degeneracy  may 
be  partially  or  wholly  lifted  by  crystal  field  effects.  Furthermore,  using  MA  to  get  Jahn- Teller 
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strengths  from  an  asymmetric  site  may  have  the  effect  of  artificially  inflating  the  MA  calculated 
values  of  the  vibronic  modes.  In  the  simulations,  the  magnitudes  of  the  vibronic  modes,  Ajr-; 
and  A 2NC,  are  consistently  larger  in  the  asymmetric  sites  than  in  the  symmetric  sites.  Typically 
the  squares  of  the  asymmetric  site  cubic  mode  parameters  were  two  orders  of  magnitude  larger, 
while  the  squares  of  the  non-cubic  mode  parameters  were  an  order  of  magnitude  larger  than 
those  for  the  symmetric  sites.  This  is  a  significant  difference,  and  overlooking  this  possibility 
may  lead  to  substantial  errors  in  analysis  of  a  system. 

The  only  possible  exception  to  the  lack  of  site  symmetry  influence  occurs  in  the  higher  order 
moments  in  the  MA  analysis.  The  ratios  of  some  of  these  are  purported  to  provide  information 
on  the  vibronic  modes  in  the  system,  as  well  as  a  second  check  on  the  magnitude  of  the  spin- 
orbit  splitting.  This  information  has  been  used  to  draw  conclusions  about  the  strength  of  the 
Jahn- Teller  active  modes  in  the  system.  There  are  a  couple  of  problems  with  this  usage.  First, 
these  values  tend  to  be  highly  random  and  second,  for  most  M/Rg  systems,  the  trapping  sites 
cannot  be  assumed  to  be  symmetric.  If  the  sites  are  asymmetric,  the  JTE  cannot  operate,  so 
the  physical  significance  of  the  higher  order  moments  is  called  into  question. 

7.  Conclusion 

When  we  initially  considered  the  experimental  results  derived  from  the  MA  method,  we  con¬ 
cluded  that  they  could  not  possibly  be  right.  This  conclusion  was  based  on  both  the  explanations 
provided  for  the  results  and  the  trapping  site  symmetry  assumptions  used  in  connection  with 
the  MA  equations.  The  theoretical  evaluations  of  the  results  either  were  not  consistent  with 
the  experimental  results  or  were  based  on  faulty  physics.  The  symmetry  assumptions  for  the 
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trapping  sites  simply  were  not  consistent  with  multiple  trapping  sites  in  a  single  rare  gas,  nor 


were  they  consistent  with  the  relative  sizes  of  the  different  metals/rare  gas  host  atoms.  However, 
MA  of  simulated  spectra  consistently  returned  the  input  SO  splitting,  providing  an  independent 
confirmation  of  the  MA  method  for  this  value  in  the  M/Rg  systems.  Further  investigation  into 
the  MA  method  showed  that  the  symmetry  assumptions  were  based  on  f-center  trapping  sites, 
and  are  not  necessary  to  MA  in  general.  Therefore,  the  simulations  and  subsequent  study  has 
shown  that  the  MA  method  is  valid  for  M/Rg  systems,  however,  the  equations  containing  the 
JT  vibrational  information  should  not  be  used. 

It  is  possible  that  the  lack  of  significant  variation  in  the  MA  generated  SO  splitting  with 
trapping  site  is  not  entirely  accurate.  Any  configuration-dependent  influences  on  the  SO  split¬ 
ting  are  not  well  determined  and  therefore  not  explicitly  included  in  the  physical  basis  for  the 
simulations.  The  importance  of  the  site  symmetry  is  more  accurately  ascribed  to  the  SO  split¬ 
ting  that  will  exist  in  the  real  system,  since  it  will  play  a  role  in  the  electron-lattice  interaction 
which  in  turn  affects  the  SO  splitting  of  the  atom  or  molecule. 

Several  tests  of  the  spectral  features  that  can  influence  the  spectra  and  therefore  the  MA 
results  were  conducted.  These  demonstrate  the  sensitivity  of  the  MA  method  to  variations  in 
the  MCD  spectra,  highlighting  the  need  for  extra  care  when  measuring  MC'D  spectra.  Further 
simulation  tests  of  the  factors  influencing  MA  results  for  M/Rg  systems  will  be  presented  in 
future  papers. 
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